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Low-temperature scanning tunneling microscopy of ring-like surface electronic 
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We report on the experimental observation by scanning tunneling microscopy at low temperature 
of ring- like features that appear around Co metal clusters deposited on a clean (110) oriented surface 
of cleaved p-type InAs crystals. These features are visible in spectroscopic images within a certain 
range of negative tunneling bias voltages due to the presence of a negative differential conductance 
in the current-voltage dependence. A theoretical model is introduced, which takes into account 
non-equilibrium effects in the small tunneling junction area. In the framework of this model the 
appearance of the ring-like features is explained in terms of interference effects between electrons 
tunneling directly and indirectly (via a Co island) between the tip and the InAs surface. 

PACS numbers: 68.35.Dv, 68.37.Ef, 73.20.At 
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After the pioneering work of Tsui [l[ , two-dimensional 
(2D) electron systems have been intensively investigated, 
in particular because of their huge application potential 
for the manufacturing of field effect transistors. Since 
then, a lot of success has been obtained with respect to 
improving the stability and the carrier mobility of the 2D 
electron gas. On the other hand, the 2D electron gas pro- 
vides an excellent playing ground for studying the phys- 
ical properties of low-dimensional systems. Up to now, 
the vast majority of experiments have been performed re- 
lying on tunneling contacts with planar geometry, where 
data is inevitably averaged and no local (on nanometer 
scale) information can be obtained. The nanometer size 
tunneling contact, which is created in scanning tunneling 
microscopy (STM) and scanning tunneling spectroscopy 
(STS) experiments, is able to directly provide such local 
information about the 2D electron gas. 



Here, we present the results of local measurements of 
the electronic properties of a 2D system, which consists of 
metal Co islands deposited in situ on the (110) oriented 
surface of InAs that is obtained by cleavage in ultra-high 
vacuum (UHV) of a single crystal. It is well known that 
by deposition of only a small amount of different adsor- 
bates onto the InAs(llO) surface, the Fermi level can be 
shifted into the conduction band by up to 600 meV (see, 
e.g., 0, H, H| and references cited therein), thus causing 
strong downward band bending and creating a 2D surface 
electron channel. 



We performed combined STM and STS measurements 
with a commercially available low-temperature STM [f| . 
The STM operates under UHV conditions at a base pres- 
sure in the 10 -12 mbar range and at temperatures down 
to 4.5 K. In order to improve measurement stability, the 
whole UHV setup is decoupled from the building by a 
specially designed vibration isolation floor. Electrochem- 
ically preformed tungsten tips [5j, which are flash an- 
nealed in situ at elevated temperatures, are used. The 
tip quality is routinely checked by acquiring atomic- 
resolution images for the "herringbone" reconstruction 
of the Au(lll) surface. 

Spectroscopic measurements were performed in the 
current imaging tunneling spectroscopy (CITS) mode @ 
with a grid size up to 200 x 200 points. Additionally, 
we relied on harmonic detection with lock-in amplifi- 
cation for "hardware" measurements of the differential 
conductance. The excitation frequency was in the range 
300 — 7000 Hz, while the excitation amplitude was in the 
range 5 — 100 mV. 

The investigated InAs samples are doped with Mn at 
a doping level of 5 x 10 17 cm -3 . At this doping level Mn 
in InAs is expected to act as a shallow acceptor with 
ionization energy around 28 meV. InAs slabs with size 
5x2x2 mm 3 were cleaved in-situ at room temperature 
in the UHV preparation chamber (base pressure is about 
5 x 10~ n mbar) that is attached to the UHV STM cham- 
ber. Co islands were deposited by means of electron- 
beam evaporation at a rate of 0.007 ± 0.001 monolayer 
(ML) per second. The evaporation rate was calibrated 
by determining with x-ray reflectometry the thickness of 
thick Co films. As discussed in more detail below, the 
evaporation rate inferred from STM images of the Co is- 
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FIG. 1: Schematic diagram of the surface band bending at the 
p-type InAs(llO) surface after Co metal cluster deposition. 
The band gap value, the width of the depletion layer and 
the position of the 2D subbands are shown. Experimental 
normalized tunneling conductance curves are shown in the 
inset both for the original, clean p-type InAs(llO) substrate 
and for the n-type Co covered substrate (Co-InAs system). 
Color online. 



lands on the InAs surface can deviate considerably from 
the calibrated rate. 

For the InAs(llO) surface, filled and empty surface 
bands are located deep inside the conduction and va- 
lence bands, respectively, and are separated by a 1.7 eV 
surface gap @, H[ for a clean surface. Under clean con- 
ditions the bands are either flat, or a weak accumulation 
layer is formed on the surface @. However, surface states 
can be strongly affected within 2 to 3 hours when the 
base pressure exceeds (2 — 3) x 10~ 10 mbar [13] • In order 
to prevent any undesirable effects of adsorbate induced 
surface degradation, we always keep the pressure below 
10~ 10 mbar in the preparation chamber of our UHV setup 
during the different sample manipulations. 

All the reported STM/STS experiments are performed 
at liquid helium temperature (~ 4.5 K). Everywhere in 
the text the tunneling bias voltage Vt refers to the sample 
voltage, while the STM tip is virtually grounded. The ex- 
periments consist of three stages. First, the clean surface 
of a cleaved p-type Mn doped InAs(llO) single crystal 
is characterized in detail by STM/STS measurements. 
Second, Co metal atoms are deposited on the InAs(llO) 
substrate that is still cold (well below room temperature) 
after removal from the STM chamber. The amount of Co 
deposited on the surface is about 0.02 ML. Depositionof 
a monolayer is according to the definition given in 
Third, the Co-InAs system is investigated by STM/STS 
measurements. 



B. 2D surface of Co-InAs 

The quality of the clean cleaved InAs(llO) surface is 
first carefully checked. Stable STM images of the surface 
reveal (lxl) relaxation with typical "butterfly" type of 
patterns appearing at the Mn doping atom sites [ll| , su- 
perimposed upon the InAs atomic lattice (see Fig. [^a)). 
Mn atoms can be visualized in the first few subsurface 
layers. Our STS measurements clearly confirm the p- 
type character of the surface conductivity (see inset in 
Fig. [TJ curve a), as expected for Mn doped InAs. Af- 
ter deposition of a very small amount of Co atoms (only 
0.02 ML) onto the surface, the general outlook of the nor- 
malized {dl / dVt) / {I /Vt) tunneling conductance curves 
drastically changes (see inset in Fig. [H curve b). The 
Fermi level Ep is now positioned about 0.12 eV above 
the conduction band bottom. Two distinct peaks are 
visible on the normalized tunneling conductance curves 
around — 0.30eV and around — 0.45eV. 

In order to provide a consistent explanation for the ob- 
served tunneling conductance curves, we need to take into 
account surface band bending (see Fig. [1]). Co clusters 
are acting as "donors" when grown on a clean InAs(llO) 
surface 3j. The Co clusters are positively charged and 
downward band bending occurs at the surface. We note 
that the average value of the positive charge residing on 
a cluster is not an integer multiple of a single electron 
charge [121 ]. 

Since our results indicate that Ep — Ec — 0.12 eV, 
with Ec the conduction band edge, and keeping in mind 
that the band gap for InAs is 0.42 eV at 4.5 K [jj], 
we then obtain a band bending at the surface that is 
as large as 0.54 eV. This gives rise to a deep quantum 
well with a few discrete eigenstates. The binding energy 
of the lowest eigenstate (obtained by solving the Poisson- 
Schrodinger equation following Q) is E ~ 0.24 eV. Con- 
sequently, the eigenstate is located almost 0.30 eV above 
Ec, or 0.18 eV above Ep. To make this state visible in 
the tunneling conductance curve, it has to be shifted in- 
side the tunneling window [Ep,Ep — eVt], and this can 
be achieved by making the tunneling bias voltage suffi- 
ciently negative. When the eigenstate Eq reaches Ep, 
a peak will appear in the normalized tunneling conduc- 
tance curve. We assume that this happens at a sample 
bias voltage of about —0.30 V for our experiments. The 
next peak appearing in the normalized tunneling conduc- 
tance curve at V f ~ —0.45 V can be assigned in a similar 
way to the fact that the eigenstate E\ falls below Ep. 
This is similar to what was described for the GaAs(llO) 
surface in [l4l ]. 

The top of valence band, Ey, is located 0.54 eV be- 
low the Fermi level. Consequently, at bias voltage Vt — 
— 0.45eV the width of the depletion layer is relatively 
small, and a tunneling current is able to flow from the 
InAs valence band to the tip. This can explain why the 
peak at V t ~ —0.45V is superimposed on a wide shoul- 
der that results from the bulk valence band states (see 
inset in Fig. [TJ curve b) 
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C. Rings around Co clusters 

According to the results of our STM measurements, 
Co atoms tend to form small clusters on the InAs(llO) 
surface despite the low substrate temperature during de- 
position. A typical STM image of a Co-InAs surface is 
presented in Fig. [DJa). From Fig. Ufa) we infer that the 
average cluster size is around 1.5 nm for a tunneling bias 
voltage of —0.65 V. It is important to explicitly refer to 
the value of the tunneling bias voltage when analyzing 
topographical STM images, because the apparent size of 
the Co clusters strongly depends on the applied bias volt- 
age. The majority of the clusters are one monolayer in 
height. A rough estimate gives an average number of 
Co atoms in a cluster (assuming a Co atom diameter of 
about 0.25 nm) that is around 30 atoms. There are 21 
Co clusters on the STM image in Fig. UK a), which has 
a size of 4.7 x 4.7 nm 2 . Consequently, the coverage that 
is inferred from Fig. [^a) is 0.068 ML, i.e., 3 to 4 times 
larger than the coverage obtained from the calibrated de- 
position rate. 

A very intriguing observation in our STM images of 
the Co islands on the InAs(llO) surface is the observation 
of ring-like features that appear around the Co clusters. 
In Fig. Ufa) one observes more or less circular depres- 
sions with sharp edges around some of the Co clusters. 
A surprising fact about these depressions is that the de- 
pressions have a different diameter at the same value of 
the bias voltage for clusters that look very similar in the 
STM image. In previous STM experiments [l5| , similar 
features were observed around metal clusters of unknown 
origin. 

In order to obtain a more detailed insight into the ap- 
pearance of the ring-like features, we performed spatially 
resolved spectroscopic measurements above the same sur- 
face area where the STM image was obtained (Fig.[^a)). 
In a first set of experiments, we relied on imaging based 
on harmonic detection with a lock-in amplifier. The re- 
sults are shown in Figs. [2Kb), (c), (d) and (f). Each image 
corresponds to a map of the tunneling current response 
at the frequency that is used for the modulation of the 
tunneling bias voltage. In a first approximation, this re- 
sponse is proportional to the value of the local density of 
states (LDOS) p according to the relation 

dI/dV t {x 7 y, V t ) oc p{x,y,E - E F = eV t ) . 

Next, we performed high-resolution CITS measurements 
with grid size of 130 x 130 points. Figure [He) shows 
a set of I(Vt) curves acquired in the CITS mode along 
the cross-section oo' indicated in Fig.[^c). Regions with 
negative slope, which corresponds to a negative differ- 
ential conductance (NDC), can be clearly observed on 
all but the two lowest curves. The latter two curves 
were acquired near point a, i.e., at points that are at 
the largest distance from the cluster. The closer to the 
cluster the I(Vt) dependence is measured, the lower in ab- 
solute value is the tunneling voltage at which the region 
with NDC occurs. This experimental finding remains 



intact for STM tip to cluster distances down to a few 
tenths of 1 nm. Comparing Fig.[^e) with Figs. [2Kb), (c), 
(d) and (f), we come to the conclusion that the distance 
d between the STM tip and the cluster determines the 
bias voltage VNDc(d) around which the NDC occurs on 
the I(Vt) curve. A dark ring with radius d, which can 
be directly related to the presence of the NDC occurring 
around V^Dcid), will appear on the dI/dVt(x, y, Vt) im- 
age taken near the bias voltage VNDc{d)- 

From our analysis of the results shown in Fig. [2] we 
conclude that (i) sharply defined dark rings appear in the 
differential conductance (LDOS) images around some of 
the clusters within a certain range of negative tunneling 
bias voltages, (ii) the size of the dark rings is shrinking 
when the absolute value of the tunneling bias voltage is 
decreased, and (iii) the LDOS has an almost constant 
value above the Co-InAs surface except for the clusters 
and the immediate vicinity of the dark rings. The only 
image, which reveals pronounced spatial variations of the 
normalized tunneling conductance, is image Fig. [2jc), 
taken at the threshold value of the tunneling bias volt- 
age Vt max that approximately corresponds to the bias 
voltage where the ring radius reaches a maximum (see 
Fig. [U(e)). We assume that the presence and the dif- 
ferent diameters of the dark rings around the different 
clusters is caused by the different bonding of the clusters 
to the InAs substrate. This assumption is supported by 
the fact that different Co clusters look quite differently 
on STM images taken at the same value of the tunneling 
bias voltage. I(Vt) curves also strongly differ when mea- 
sured above different clusters. Curves, which are mea- 
sured above clusters that are surrounded by a ring-like 
feature, reveal the presence of a narrow plateau (width of 
about 0.1 V) around the Fermi level. On the other hand, 
clusters, which are not surrounded by a ring-like feature, 
reveal the presence of a wider plateau in the I(Vt) curve 
(width is 0.5 — 0.7 V). These findings clearly require fur- 
ther investigation. 

The maximum observable ring radius is determined by 
the fact that the region with negative differential con- 
ductance (see Fig. [2Je)) becomes masked by the rapidly 
growing tail of the I(Vt) characteristics at more negative 
bias voltages. For the Co clusters in the selected area 
this masking occurs for voltages V t < VtuAX — —0.7 V 
while for other areas the value of VtMAX will be different. 
In this bias voltage range the behavior of the tunneling 
conductance becomes much more complicated. As can 
be seen in Fig. [2Kb) , ring- like protrusions are observed in 
the LDOS image rather than ring-like depressions. Some 
of the rings become partly dark and partly bright. For 
these more negative bias voltages, the exact relation be- 
tween the local slope of the I(V t ) curve and the depres- 
sion of the tunneling current becomes important. More 
subtle details, such as the presence of neighboring clus- 
ters, impurity atoms or surface defects, will then have a 
pronounced influence on the imaging of the clusters. 

The large grid size that has been used for acquiring the 
CITS data allows us to present the data in another inter- 




FIG. 2: (a) Experimental constant current topographical STM image of the Co-InAs(llO) surface with a Co coverage of about 
0.02 ML obtained at Vt = —0.65 V and / = 400 pA. The image size is 47 x 47 nm 2 . (b), (c), (d) and (f) are maps of the measured 
differential conductance dl /dVt(x,y, Vt) for the same surface region as in (a) at different tunneling bias voltage. The tunneling 
current setpoint is It = 150 pA, while the tunneling bias voltage is indicated at the bottom of the maps. Ring-like features, 
whose appearance is strongly affected by the applied negative bias voltage, can be observed, (e) Set of I(Vt) curves extracted 
from the CITS data (see text) along the line aa' (see (c)). The different curves are offset for clarity. Regions with negative 
differential conductance are clearly visible. The value of the tunneling bias voltage corresponding to the maximum observable 
ring radius (see text) is marked as VtMAX- The ring in (c) that is below the line aa' has almost reached its maximum radius. 



esting way. Figure [J2a) gives a 2D map of the variation 
of the normalized tunneling conductance data along the 
cross-section bb' in Fig. [^a). The horizontal axis corre- 
sponds to the spatial coordinate, while the vertical axis 
corresponds to the tunneling bias voltage, and the color 
intensity corresponds to the measured value of the nor- 
malized tunneling conductance (LDOS). In Fig.^b) we 
show (dl I 'dVtj I '{I /Vt) at the three points that are marked 
by the crosses on the cross-section bb' in Fig. HJa). These 
three points (1, 2, and 3) are marked in Fig. 03a) as well. 

From Fig. [3ja) it is clear that a Co cluster affects the 
LDOS in the conduction and valence bands only very lo- 
cally. Except for a limited voltage range near Vt = 0, 
the influence of a cluster on the LDOS is very rapidly 
decaying outside the cluster. Inside the cluster, the con- 
ductivity in the valence band is considerably suppressed. 
Some of the energy levels of the Co cluster become visible 
in the voltage range corresponding to the band gap and 
the conduction band. In Fig. [3^a) there are two tilted 
dark lines that approximately start from the peak in the 
normalized tunneling conductance of the cluster and that 
extend over a distance of about 8 nm into the defect free 



surface area. These two dark lines directly reflect the 
presence of the dark rings around the clusters as well as 
the dependence of the ring diameter on the tunneling bias 
voltage, which is clearly non-linear. 

The two peaks in the normalized tunneling conduc- 
tance at —0.3 V and at —0.45 V (see curve (a) in the inset 
to Fig. [J) also appear in Fig.[3^a) in the defect free surface 
area. The peaks are visible as perfectly straight, horizon- 
tal stripes. This implies that there occurs no noticeable 
change in the "rigid" (both bands move together) band 
shifts relative to the original bands of the p-type InAs. 
This is rather surprising, because a charged cluster and a 
charged impurity atom, separated by a distance of about 
20 nm, are surrounding the scanned area, and one might 
therefore expect changes to occur in the measured band 
shift. The absence of changes in the rigid band shift may 
be accounted for in the following way. First, both the 
Co cluster and the Mn impurity atom can be treated 
as point charges. It is known that a point charge po- 
tential does not strongly affect the band bending, but 
rather gives rise to fine structure in the LDOS [Tj|]. Sec- 
ond, the variations of the LDOS around the Co cluster 
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are strongly confined in space, while the Mn impurity 
atom alters the LDOS mostly near the top of the valence 
band. Hence, the conduction band remains almost undis- 
turbed, with the conduction band bottom Ec appearing 
as a straight line in Fig. [Ha). When the tunneling bias 
voltage Vt , which is sweeping during the measurements of 
the tunneling conductance, reaches the band gap region 
and goes down to the top of the valence band, the energy 
levels E , Ei of the 2D quantum well become populated. 
Taking into account the very high electron mobility along 
the surface [jq |. we may assume that in this bias range 
perfect screening by the 2D electron gas is taking place, 
hiding other details in the electronic structure. 

On the other hand, when the 2D subband E\ is shifted 
below the Fermi level, the tunneling bias dependent band 
bending noticeably slows down. This can be accounted 
for by the overlapping of the higher 2D subbands, which 
form a relatively wide continuum of states, thus pinning 
the Fermi level in this continuum. In this case the sur- 
face band structure "stabilizes", i.e., an increase of the 
absolute value of the negative bias voltage leads only to 
slight changes in 2D subband energy. 

Another important observation related to Fig. OJa) is 
the influence of the Mn impurity atoms on (below) the 
Co-InAs(llO) surface on the tunneling spectra. In the 
vicinity of a Mn atom the LDOS is affected deep in- 
side the conduction and valence bands (not visible in 
Fig. El^a)). Additionally, an intense peak appears in 
Fig. GD(a) near the valence band top and most probably 
this peak reflects the energy position of the Mn acceptor 
band. We note that the spatial extent of this peak is 
relatively large, about 7nm. 



III. DISCUSSION 

In order to explain the observed effects consistently 
one needs to take into account the presence of non- 
equilibrium and interference effects in the small tunneling 
junction [17], which can be described using the Keldysh 
diagram technique. 

We need to explain two main peculiarities that are re- 
lated to the observation of the ring-like features. First, 
the presence of the NDC region in the tunneling conduc- 
tivity curves has to be accounted for. Second, the limited 
spatial extent of the circular features requires an expla- 
nation. Accordingly, this discussion section has been di- 
vided into two subsections. 
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FIG. 3: (a) 2D representation of the normalized differential 
tunneling conductance, extracted from the CITS data (see 
text) along the line bb' indicated in Fig. [2ja). The horizontal 
axis corresponds to the spatial coordinate, while the verti- 
cal axis corresponds to the tunneling bias voltage. The color 
contrast in the image corresponds to the value of the differ- 
ential tunneling conductance as given by the color scale on 
the right hand side, (b) Three normalized differential con- 
ductance curves measured at the three points (1, 2, and 3, 
see (a)) that are marked by the crosses along the line bb' in 
Fig-EJa). The curves are measured above a cluster (1), above 
a defect free surface area (2) and above a Mn impurity atom 
(3). 



where 

Ho = ^2(s q " M)Cq CT c qCT + ^2(e p -n- eV)c+ a c pa 

q,cr p,o" 

describes the non-interacting electrons in the two elec- 
trodes: tip states (per), and bulk states (qer). The part 
H c i corresponds to the electron states in the cluster: 



H c i = e a ^2 ' 



A. Negative differential conductance 



We start from the following model Hamiltonian H: 



H — H + H tun + H c i + Hqw 



where the operator a a (a+) destroys (creates) an elec- 
tron in the cluster with spin a. The part Hqw describes 
the electron states in the quantum well resulting from 
the Co cluster deposition and the STM tip induced band 
bending: 

Hqw = £kmo-Ckmo-c£ TO(7 , 

k,m,cr 
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FIG. 4: Schematic representation of the tunneling processes 
occurring in the system consisting of the STM tip, a Co clus- 
ter and the InAs(llO) substrate. e a is the Co cluster energy 
level that is participating in the tunneling. Eq is the first 
2D subband, which has a non-zero but narrow width, as in- 
dicated by the gray stripe. A gray stripe is also used to mark 
the presence of the higher 2D subbands, which form a contin- 
uum of states due to their broadening. The system is shown 
at non-zero bias voltage when e a and Eo are aligned and a 
dark ring appears around the Co cluster in the maps of the 
differential tunneling conductance. 



where c^ mcr (c^ ma .) destroys (creates) an electron with 
spin a in the quantum well, m is the number labeling the 
2D bands and k is the 2D wavevector. H tun is responsible 
for tunneling transitions from the quantum well and from 
the cluster into each of the electrodes (tip or substrate): 



Hf, 



E 



p.cr 

p,k,m,a" 



q,k,m,cr 



T pa , Tk m , T q and W are the relevant tunneling matrix 
elements: Tk m for tunneling from the cluster level to the 
quantum well subbands, T pa for tunneling from the clus- 
ter level to the STM tip continuum, T q for tunneling from 
the quantum well level to the bulk of the semiconductor, 
W for tunneling from the quantum well level to the STM 
tip. 

The differential tunneling conductance can then be 
written as [13, LL9( 



dI r , 2 V0V7feTVfe(l - IV fe ) 2 R 
-ocr, fe + ^—^ ™^eG (eV) 



2 7p7fc((l + Tv k f - ATv k cos 2 ip) 



lmG R {eV) 



(7 P + 7fc)(i + r^ fe ) 3 

2r^ fc (7fc +7 P ) 

(i + iv fe ) 3 



ImG R (eV) , (1) 



were Uk, Vp and v q are the density of states in the quan- 
tum well, the STM tip and the bulk semiconductor, re- 



spectively. The cluster density of states is proportional 
to Im G R (ui). 

The different tunneling rates in Eq. ([1]) are defined 
as follows (see Fig. [3J: Jk = T^ m v km is the tunneling 
transition rate from a cluster level to the quantum well 
subbands, 7 P = T 2 v p is the tunneling transition rate from 
a cluster level to the STM tip, -f q = T^v q is the tunneling 
transition rate from the quantum well to the bulk of the 
semiconductor, and 7 = W 2 v p is the tunneling transition 
rate from the quantum well to the STM tip. Finally, T 
is defined as T = 77,7/(7 + 7 9 ). We omit the subband 
index to, because only the lowest subband plays the most 
important role in our model. 

As stated above, a Co cluster acts as a donor, and 
consequently it is positively charged. In case of non- 
equilibrium, the energy of a cluster level s a , which is 
participating in the tunneling (see Fig. [4j, depends on 
the electron filling number < n a > of this level: 

£ a = £q + U(n - < n a >) , 

where eo is the unperturbed position of the cluster energy 
level, no is the equilibrium cluster level filling number, 
and U is the Coulomb interaction energy of the localized 
charges (l7| . 

For STM junctions we can assume that Tvk « 1, 
because V is of order of tunneling rates and Vk ~ 1/-D, 
where D is the width of conduction (valence) band in 
In As. Equation |T]) then reduces to 



A(^) oc 2 I Im G R {u) \ l^p^k+Jhl 



7p + Ik 



cos if , (2) 



where A refers the the extra contribution to the tunneling 
conductance resulting from the non-equilibrium effects. 
Here tp is the pase angle, ascuired due to difference of 
wavefunctions of direct (STM tip-2D subband) and in- 
direct resonant tunneling. The complicated coefficient 
in front of Im G r (lo) and the presence of Re G R (ui) 
in Eq. @ also results from interference between direct 
tuneling from the STM tip to a quantum well subband, 
and indirect resonant tunneling via the cluster level. 
Consequently, the involved tunneling rates determine the 
sign and the magnitude of the additional current result- 
ing from the non-equilibrium and interference processes. 

For each tip position the value of 7 P is fixed, while 7^ 
can be varied by changing the applied bias voltage. 7& 
depends on the relative position of the cluster level s a 
and of the first 2D subband Eq. Both energies are bias 
dependent: e a trough < n a > and Eq through the band 
bending. It is important to note that the 2D subband 
Eq is relatively narrow. Electrons can be injected into 
this subband either from the STM tip (7) or from the 
Co cluster (7/.). In both cases the wavevector compo- 
nent parallel to the surface, fen, is small. In close analogy 
with angle resolved photoemission spectroscopy experi- 
ments [20|, only a narrow range of energies is available 
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for tunneling into the 2D subband Eo (see Fig. [4]). In gen- 
eral, the fc| I component needs to be taken into account in 
STM experiments 

We assume that the unperturbed Co cluster energy 
level Eo (at zero bias voltage) is positioned close to Ep 
below the first 2D subband Eo, i.e, the two energy levels 
are not aligned. This assumption is strongly supported 
by the experimental STS data (see Fig. [3] for details). 

We now proceed with a more detailed analysis of the 
bias voltage dependence of A(dI/dV). The first term in 
Eq. ([2]), which varies proportional to the cluster density 
of states, can change its sign. If jpjk « r^ fc (7 p + j k ) 2 
(let us assume j p << jk) the sign of the coefficient that 
appears in front of Im G r (uj) is negative. For the as- 
sumed conditions, the coefficient that appears in front 
of Re G r (uj) is much smaller than the coefficient next 
to Im G R (ui). The tunneling conductance curve is an 
inverted Lorentzian-shaped curve that is slightly asym- 
metric because of the Re G r (oj) contribution. The posi- 
tion of the dip in the tunneling conductance curve can be 
connected to the e Q (< n a >) dependence. The closer to 
the cluster the STM tip is, the higher is the value of j p 
and the lower is the value of the non-equilibrium filling 
< n a >. The cluster energy level e a is pushed higher, 
and the dip in the tunneling conductance curve can be 
observed at smaller absolute values of the applied bias 
voltage (see Figf2|). The dark rings that appear in the 
STS images shrink when the absolute value of the nega- 
tive tunneling bias decreases within a cluster dependent 
range of negative voltages. 



B. Spatial extent of ring-like features 

The most probable physical mechanism causing the 
formation of narrow (of the order of the interatomic dis- 
tance) ring-like features in real space is some kind of res- 
onant tunneling. When the STM tip changes its spatial 
position, localized energy levels that participate in the 
tunneling processes are coming out of resonance or even 
are falling outside the tunneling window [Ep, Ep — eVt]. 
We find that the simple cluster charging model used in 
[l5L |22| is not adequate for describing our Co-InAs sys- 
tem. This is mainly due to the fact that on-site Coulomb 
repulsion for the cluster is the strongest bias dependent 
effect. This leads to variations of the cluster level filling 
number, and consequently to changes in the cluster level 
energy. Starting from this finding we are able to ana- 
lyze in the framework of our model how spatially narrow 
circular features can be formed around the Co clusters. 

For each value of the tunneling bias voltage Vt the 
value of 7 P can be altered by changing the STM tip to 
cluster distance. Consequently, < n a >, which depends 
on the ratio "fk/"fp also changes, giving rise to a shift of 
£a{< n a >) and to variations of 7&. When moving the 
STM tip farther away from the cluster, j p is decreasing, 
leading to an increase of < n a > and to a lowering of 
e a . At a certain distance d max , the cluster energy level 



£ a {< n a >) starts to becoming out of resonance with 
the narrow 2D subband Eo- Since there is now an align- 
ment with the tail of the broadened Eo level, the value 
of the tunneling current drops. At the same time, the 
value of 7^ decreases, resulting in a decrease of the non- 
equilibrium occupation < n a >, and consequently e a is 
pushed higher. The cluster level e a then jumps to its 
equilibrium position So- We assume that at the STM tip 
to cluster distance d max the tunneling rate j p , which ex- 
ponentially decreases with increasing distance, becomes 
negligible (7 p 1=3 0). The STM tunneling system does 
no longer feel the STM tip - cluster interaction, and the 
value of the tunneling current returns to its "normal" 
value. In other words, when at a certain tunneling bias 
voltage in the range where a dark ring appears in the 
LDOS around the Co cluster, the STM tip to cluster 
distance d reaches the value d max , some kind of positive 
feedback, caused by mutually coupled variations of 7 p , 7^. 
and < n a > , switches on. This positive feedback abruptly 
shifts the cluster energy level out of resonance with the 
2D subband. This can account for the very small spatial 
width of the rings that appear in the LDOS as well as 
for the non-linear dependence of the ring diameter on the 
STM tip to cluster distance (see Fig. [3]). 

We note that for specific tunneling conditions the first 
term in Eq. ^ can become positive, causing the obser- 
vation of bright rings or protrusions in the images of the 
differential conductance. A dark ring can turn into a 
bright one at high values of the tunneling bias voltage, 
as can be seen in Fig. (Hb). 

Our theoretical model does not include any assump- 
tions concerning the physical origin of the localized states 
through which resonant tunneling can occur. Therefore, 
this model should also be applicable to the case of an 
individual impurity atom near the surface. Ring-like fea- 
tures around individual impurities were indeed observed 
in our experiments (to be published) as well as on the 
GaAs(llO) surface [j^. 



IV. CONCLUSION 

In conclusion, relying on combined STM and STS mea- 
surements at low temperature, we have been able to iden- 
tify the presence of ring-like features around Co metal 
clusters on a p-type InAs(llO) surface that is prepared 
by in situ cleavage of a single crystal in UHV. The ring- 
like features become clearly visible in the maps of the 
differential normalized conductance for a certain range 
of negative tunneling bias voltages, due to the presence 
of a region with negative differential conductance in the 
I{Vt) dependence. The diameter of the rings is decreas- 
ing when decreasing the absolute value of the tunneling 
bias voltage. The possibility to observe the ring-like fea- 
tures around a Co cluster depends on the bonding of the 
cluster to the InAs(llO) surface. 

A theoretical model was developed, which takes into 
account the non-equilibrium effects that occur in the 
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small STM tunneling junction. The model we propose 
allows us to explain the main experimental findings. In 
the framework of the model the appearance of the ring- 
like features structures is accounted for by the interfer- 
ence of direct (between the STM tip and the QW) and 
indirect tunneling via a Co cluster. 
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